Cohesion of BaReHg and BaMnH 9 : Density Functional Calculations and Prediction of 

(MnHg) 2 " Salts 
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Density functional calculations are used to calculate the structural and electronic properties of 
BaReHg and to analyze the bonding in this compound. The high coordination in BaReHg is due 
to bonding between Re 5d states and states of d-like symmetry formed from combinations of H s 
orbitals in the Hg cage. This explains the structure of the material, its short bond lengths and other 
physical properties, such as the high band gap. We compare with results for hypothetical BaMnHg, 
which we find to have similar bonding and cohesion to the Re compound. This suggests that it may 
be possible to synthesize (MnHg) 2- salts. Depending on the particular cation, such salts may have 
exceptionally high hydrogen contents, in excess of 10 weight %. 

PACS numbers: 71.20.Lp,71.20.Be,61.50.Lt 
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I. INTRODUCTION 

BaReHg is the prototypical member of a family of 
hydrides, based on (ReHg) 2 ~ and (TcHg) 2 ~ structural 
un j£ S| i i 2 4 3 j 4 i 5 4 6 These compounds are of interest from 
a fundamental point of view* because of the unusual 
coordination of the transition metal atoms, which are 
surrounded by 9 hydrogen atoms, with relatively short 
metal-H bond lengths and are in an high formal valence 
state of 7. These compounds are also of practical inter- 
est because of the high hydrogen to metal ratio of 4.5 in 
BaReHg . Although no Mn based examples of these com- 
pounds have been synthesized to date, the hypothetical 
3d analogue, MgMnHg, would have a hydrogen content 
in excess of 10 weight percent. The purpose of this pa- 
per is to analyze the electronic structure of BaReHg in 
order to understand its bonding and the prospects for 
synthesis of Mn based analogues. This follows a previous 
electronic structure calculation by Orgaz and Gupta, 8 
done using the X-ray crystal structure of Ref. and the 
Xa method. The present calculations were done using a 
general potential self-consistent method, with calculated 
H positions, and allow us to present a more accurate elec- 
tronic structure and detailed analysis of the bonding. 

The crystal structure (spacegroup P63/mmc), which 
was determined using X-ray diffraction by Stetson and 
co-workers^ is depicted in Fig. ^ It consists of alternat- 
ing triangular layers of Ba and Re stacked along the o 
axis, such that the Re are at the center of trigonal prisms 
formed by the Ba. Each Re is coordinated by three H 
atoms in its own plane (the HI) positions, and six other 
H atoms (the H2) above and below the plane. These H2 
atoms form a trigonal prism. Thus the H atoms are ar- 



ranged in tricapped trigonal prisms around the Re atoms. 
The unit cell contains two (ReHg) 2- units, stacked so 
that the orientation of these units alternates along the 
hexagonal c-axis. 



II. METHOD 

The calculations reported here were done within 
the local density approximation (LDA) using the gen- 
eral potential linearized augmented planewave (LAPW) 
method. 9 Local orbitals^ were used for the semicore 
states (Re 5s and 5p, and Ba 5s and bp), and to relax 
the linearization of the d bands. LAPW sphere radii of 
2.0 ao and 1.1 ao were used for the metal and H atoms, 
respectively. We used well converged basis sets consist- 
ing of more than 3200 LAPW functions and local orbitals 
for the two formula unit primitive cell. We tested larger 
basis sets, but found no significant changes in the re- 
sults. We also tested various Brillouin zone samplings 
for the iteration to self-consistency, but found that the 
results were already converged at a sampling of 6 special 
k-points in the irreducible wedge of the hexagonal zone. 
No doubt this reflects the large band gap insulating char- 
acter of the compound. The calculations for hypothetical 
BaMnHg were done with similar parameters, except that 
the Mn and H radii were 1.8 ao and 1.05 ao, respectively. 



III. STRUCTURE 

While the lattice parameters, metal positions, and ba- 
sic structure determined by Stetson and co-workers (Ref. 




FIG. 1: (color online) Structure of BaReHg, showing Ba and 
ReHg layers stacked along the c-axis (top) and the structure 
of the hexagonal planes (bottom). The HI and H2 sites are 
shown by small light and dark red spheres, respectively, Ba by 
large light green spheres, and Re by large dark blue spheres. 
The atomic coordinates are those obtained by LDA structure 
relaxation. 



are undoubtedly correct, it is difficult to accurately de- 
termine H positions in a complex compound containing 
heavy atoms without neutron scattering. Accordingly, 
we used the experimental lattice parameters (et=5.287A, 
c=9.323A) here, but determined the internal coordinates 
of the H atoms by energy minimization within the LDA. 
The calculated structural parameters are given in Ta- 
ble U along with those of hypothetical BaMnHg, which is 
discussed later. As may be seen, the tricapped trigonal 
prisms are quite compact and regular in that the Re-Hl 
and Re-H2 bond lengths are short and similar, although 
there is a compression along the c-axis, as may be seen 
from the fact the the H2-H2 neighbor distance on the 
top of the trigonal prism is longer than the H1-H2 dis- 
tance between the prism and cap hydrogen sites. It is 
also notable that the H1-H2 distance of 1.92A is very 
short compared with most metal hydrides, which gener- 
ally have H - H distances larger than 2.1 Aiii*i£ These 
short H-H distances within the (ReHg) 2- units suggest 
that direct H-H interactions may be important in forming 
the electronic structure, which in fact is what we find. 



FIG. 2: Calculated band structure of hexagonal BaReHg, us- 
ing the relaxed crystal structure. 
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FIG. 3: (color online) Calculated electronic density of states 
(DOS) of BaReHg on a per formula unit basis. The top panel 
shows the total, Re d and Ba contributions, as defined by 
projections onto the LAPW spheres. The bottom panel shows 
the projections of HI s and H2 a character. Note that for 1.1 
ao H LAPW spheres approximately 2/3 of the H a charge for 
a H atom would lie outside the sphere. 
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TABLE I: Calculated structural parameters of PQ^/mmc 
BaReHg and hypothetical BaMnHg, obtained within the 
LDA. The lattice parameters were fixed at the values 
a=5.287A, and c=9.323A, as reported in Ref. for BaReHg. 
For hypothetical BaMnHg we used a=5.067A, and c=8.883A 
(see text). The coordinates of the atoms are Ba (2a): (0,0,0), 
Re/Mn (2c): (1/3,2/3,1/4), HI (6ft,): (sci,2a;i,l/4) and H2 
(12k): (x2,2x2,Z2). Distances, d, are the shortest bonds of a 
given type. 



parameter 


BaReHg 


BaMnHg 


Xl 


0.151 


0.161 


X2 


0.466 


0.460 


Z2 


0.122 


0.129 


d(Re/Mn-Hl) 


1.67A 


1.51A 


d(Re/Mn-H2) 


1.71A 


1.54 A 


d(Hl-H2) 


1.92A 


1.73 A 


d(H2-H2) 


2.10A 


1.92 A 


d(Hl-Hl)* 


2.39A 


2.44A 


d(Ba-Hl) 


2.71A 


2. 63 A 


d(Ba-H2) 


2. 89 A 


2. 80 A 


* between different (ReHg) 2 units. 



IV. ELECTRONIC STRUCTURE 

Our main electronic structure results are given in Figs. 
121-01 The calculated band structure and electronic den- 
sity of states (DOS) are shown in Figs. [5] and [21 respec- 
tively. Fig. |3]also shows projections of the DOS onto the 
LAPW spheres. Electron counts corresponding to inte- 
gration of the DOS and projections are shown in Fig. 0] 
As in the calculation of Ref. la, we obtain a wide band gap 
insulator, but the bands and the structures and positions 
of the features in the DOS are significantly different. We 
obtain a direct band gap of 3.58 eV at T, which may 
be an underestimate, as is common in LDA calculations. 
In any case, this large band gap is consistent with the 
observed transparent nature of the material, and also in- 
dicates the nature of the band formation. In particular, 
it is an extremely large value for a crystal field gap, and 
also is not consistent with an ionic gap, since in a scenario 
with H" and Re 7+ ions, one would expect low lying Re 
d levels and therefore a small gap due to the high Re va- 
lence. Such an ionic scenario is also unlikely considering 
the short H-H bond lengths in the (ReHg) 2 ~ units. 

Examining the band structure in greater detail, we 
note that there are three manifolds of states. These are 
a narrow set of two bands (one band per formula unit) 
extending from -6.2 eV to -5 eV (relative to the valence 
band maximum), a manifold of 16 bands extending from 
-3.3 eV to the valence band maximum and a broad set 
of conduction bands. While valence bands are often less 
dispersive than conduction bands in semiconductors, this 
electron-hole asymmetry is particularly noticeable in the 
BaReHg. As may be seen from the projections, of the 
DOS there is no appreciable Ba contribution to the va- 
lence bands, and therefore, as expected, Ba occurs as 
Ba 2+ . There is, however, Ba d character in the upper 
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FIG. 4: (color online) Integrated DOS of BaReHg on a per 
formula unit basis. The top panel shows the total, Re d and 
Ba contributions, as in Fig. [3] The bottom panel shows the 
projections of HI s and H2 s character. The integrations are 
on a per formula unit basis. 



conduction bands starting at ~ 7 eV, and the Ba / res- 
onance is visible from ~ 10 - 11 eV. 

With the exception of the lowest, split-off peak in the 
DOS, which is H s in character, the remainder of the DOS 
shows a mixture of Re d and H s character, demonstrat- 
ing covalency between Re and H. This is distinct from 
what is seen in NaAlH4, for example. 13 

The projected densities of states and electron counts 
are based on integration within the LAPW spheres. In 
the LAPW method these spheres are constrained to be 
non-overlapping, and in addition they must be chosen so 
that core states are contained within them. These con- 
straints require the use of rather small spheres at the H 
sites. Here we used rn = 1.1 ao. This is significant for 
the interpretation of the projected DOS, because most 
of the H s charge lies outside this radius. In particular, 
for a neutral non-spin-polarized H atom in the LDA only 
0.329 electrons are contained within 1.1 ao, and 0.305 
electrons within 1.05 ao. This factor needs to be kept 
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in mind when analyzing the DOS. The advantage of the 
small H LAPW sphere is that the H wavefunction near 
the nucleus is much less sensitive to the environment, 
which means that charge within this size sphere is al- 
most entirely due to the H s states, with very little con- 
tamination by charge entering from neighboring atoms. 
We verified this by integrating the p projection in the H 
spheres, which would come from tails of neighbors, and 
found less than 0.02 electrons, supporting this conclu- 
sion. Thus, integration of the H s projection can provide 
a useful measure of the H valence, when normalized by 
the charge of a neutral H. The charge obtained in this 
way corresponds to H 05_ for both the HI and H2 sites. 
This analysis cannot be easily applied to the conduction 
bands because there is no plateau at the top of the Re d 
derived DOS, presumably because of contributions from 
the free electron like states, which would provide H 2s 
and other contributions. Nonetheless, by this measure, 2 
electrons per H would be included by integrating up to 9 
eV, while integrating to 10 eV, where the Re d character 
becomes small, yields 2.1 electrons. In any case, these 
numbers indicate substantial covalency between Re and 
H. This is also consistent with the large bandgap, which 
comes from the bonding - antibonding splitting associ- 
ated with the covalency. 



V. BONDING 

To understand the bonding, we examined the symme- 
try of the H s states when expanded about the Re site, 
using the projections of the H tails into the Re sphere. 
Based on this, the lowest split-off peak containing 2 elec- 
trons comes from the symmetric combination of the H 
s states, and is s-like at the Re site. The second peak, 
which contains 16 electrons per (ReHg) 2 ~ unit has both 
p-like and d-like combinations of the H s states. Con- 
sidering the short H-H distance in the Hg cage, and the 
simple electronic structure of H, the electronic states of 
a H shell should be ordered according to the number of 
nodes, yielding, in order, a fully symmetric state (s-like), 
three states with one node (p-like) and five states with 
two nodes (d-like). In fact, calculations for free Beg and 
Mgg clusters which have 18 electrons and single s or- 
bital ions, show that these clusters are exceptionally sta- 
ble and have a tricapped trigonal prismatic ground state 
structure, similar to 18 electron (ReHg) 2 ~. Furthermore, 
the electronic structure^ shows clearly separated groups 
of states, which are, in order, s-like, p-like and d-like. 
The energy splittings of the states within the groups due 
to the non-spherical, D^h, symmetry are small compared 
to the splittings between the groups. This strong shell 
structure of the orbitals is associated with the stability 
of these clusters and their shape, as discussed by Reim- 
mann and co-workers^£ and may be expected to be a 
general tendency in 18 electron, 9 site clusters with only 
s-orbitals. 

This forms the basis for understanding the electronic 
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FIG. 5: Schematic depiction of the formation of the BaReHg 
electronic structure from the levels of a tricapped trigonal 
prismatic Hg cluster (see Ref. HEF) and the Re d orbitals. The 
levels of the cluster are labeled according using Dzh represen- 
tations, and also s, p and d according to their approximate 
angular momentum character when expanded around the cen- 
ter of the cluster. Note that all the H states shown are derived 
from different linear combinations of H Is orbitals. 



structure and bonding of BaReHg, as schematically il- 
lustrated in Fig. Specifically, the d-like states of the 
Hg shell have angular character about the Re site that 
favors bonding with the Re d-states. This yields five 
bonding states, which are a bonding combinations in the 
radial direction, and five antibonding states, which have 
nodes between the Re and the H shell. This provides an 
explanation for the asymmetry of the valence and con- 
duction bands, since the antibonding combinations will 
be directed away from the (ReHg) 2 " clusters, leading to 
greater dispersion than for the bonding combinations. 

The strong bonding of the c?-like combinations of H s 
states on the Hg shell and Re d states leads to a large 
band gap, and shifts the bonding states so that they 
are in the same energy region as the p-like combina- 
tions of shell states, which are non-bonding with Re d 
orbitals. Thus, the 16 electron main peak in the valence 
DOS comes from 10 electrons in bonding Re d - shell d 
states, and 6 electrons in shell p-like states. Consider- 
ing the large bonding - antibonding splitting of ~ 8 eV 
as measured by the energy separation of the centers of 
the valence and conduction Re d DOS, there should be 
approximately equal mixtures of Re d and H s character 
in both the bonding and antibonding states. This would 
yield an electron count of 13 H Is electrons in the va- 
lence bands (including the s-like and p-like shell states) 
and therefore a hydrogen valence of H 044 ~ , in reasonable 
accord with the estimate obtained from integrating the 
projected DOS. 
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VI. HYPOTHETICAL MANGANESE BASED 
COMPOUNDS 

The bonding mechanism discussed above is quite gen- 
eral in principle, depending on a proper electron count 
and the formation of strong bonds between d orbitals 
and appropriate combinations of H s orbitals. The d or- 
bitals of 3d elements are much less extended than those 
of Ad and bd elements. Thus, as is well known, 3d com- 
pounds typically have smaller band widths and reduced 
covalency with ligands than 4d and 5d analogues. This 
might be used to provide an explanation for why salts 
containing (MnHg) 2- have not been reported. Nonethe- 
less, considering that the bonding in BaReHg appears to 
be very strong, based on the large bonding - antibond- 
ing splitting and the fact that the phase competes with 
the very stable hydride BaH 2 , it is of interest to study 
the Mn based analogue. We note that the binding en- 
ergy of MgBi2 is less than half of that of Bafi2 , and that 
the binding energy of BeBi2 is even lower, and so even 
with reduced stability of the anionic cluster, salts may 
be formed with appropriate cations. 

In order to investigate the electronic structure of hy- 
pothetical BaMnHg, we reduced the a and c lattice pa- 
rameters of the Re compound based on the difference 
in the covalent radii of Re and Mn (O.llA), to obtain 
a=5.067A, and c=8.883A. This corresponds to a volume 
reduction of 12.5% relative to BaReHg. Using these lat- 
tice parameters, we relaxed the H positions in the unit 
cell. The resulting H coordinates and bond lengths are 
give in Table [I] Remarkably, the relaxed structure main- 
tains a cage of nine H atoms at approximately the same 
distance from Mn, and with short H-H distances. 

As may be seen from the DOS and projections (Fig. 
El, the electronic structure is qualitatively similar to 
BaReHg, although the detailed peak structure of the 
main valence band DOS differs and the split-off H state 
moves to lower energy. This larger splitting between 
the shell s-like state and the shell p-like states is what 
is expected based on the smaller H-H distances in the 
more compact Hg cage of the Mn compound, since this 
will yield larger hopping matrix elements and therefore 
stronger shell structure. As in the Re compound, integra- 
tion of the H s projection of the DOS, normalized by the 
amount of charge within a 1.05 ao sphere for a neutral 
non-spin-polarized H atom in the LDA, yields a hydrogen 
valence, H 0,5_ . Furthermore, although the band gap of 
BaMnHg is reduced to 3.0 eV, the bonding - antibonding 
splitting as measured by the centers of the Mn projection 
of the DOS is nearly the same as in the Re compound. 

The implication of the above results is that salts con- 
taining (MnHg) 2 ~ may be realizable. Synthesis may, 
however, be difficult because of competition from sta- 
ble hydride phases. As noted, MgH 2 is less stable than 
BaH2. However, MgH2 still has a binding energy of 75 
KJ/mol H2, and furthermore Mg forms other stable hy- 
drides with Mn, such as Mg 3 MnH 7 iLi^ 

The LDA often overbinds solids. Nonetheless, it is of 
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FIG. 6: (color online) Calculated electronic density of states 
(DOS) of hypothetical BaMnHg on a per formula unit ba- 
sis, as m Fig. 03 Note that the H LAPW sph ere radius for 
BaMnHg was 1.05 ao. 

interest to compare the heats of formation of BaReHg 
and hypothetical BaMnHg. We calculated these values 
by subtracting calculated total energies of the elemental 
metals and the H2 molecule from the total energy of the 
compounds. For Mn, we used the fee structure and made 
a correction of 5 mRy per atom, which comes from en- 
ergy difference between the fee structure and the ground 
state non-collinear antiferromagnetic a-Mn, as obtained 
by Hobbs and Hafner^S We made no correction for zero 
point energy. 

With the LDA value of the energy of H2, as obtained 
in a supercell calculation, E(H.2) = -2.294 Ry, the result- 
ing formation energies are -99 KJ/mol H2 for BaReHg 
and -86 KJ/mol H2 for hypothetical BaMnHg. The main 
uncertainties are (1) the crystal structure of BaMnHg, 
which if different from that assumed would lead to a 
higher stability, (2) the neglected zero point motion, 
which typically affects the cohesive energies by less than 
20 K J /mol H2 and here would likely destabilize the com- 
pounds as the H- metal and H-H bonds are short, suggest- 
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ing stiff phonons, and (3) the use of the LDA, which based 
on calculations for other metal hydrides jSLSLSSi overesti- 
mates stability typically by ~ 20 KJ/ mol H2, though 
this varies from compound to compound. In any case, 
based on these results, we conclude that BaReHg and 
hypothetical BaMnHg have similar stabilities, assuming 
that the decomposition is to similar products, e.g. Bati2 
and Re/Mn metal. 

VII. SUMMARY AND CONCLUSIONS 

Electronic structure calculations show that the high 
coordination in BaReHg is due to bonding between Re 
5d states and d-like states formed from combinations of 
H s orbitals in the Hg cage. This explains the structure 
of the material, its short bond lengths and other physical 
properties, such as the high band gap. Similar bond- 
ing is found in hypothetical BaMnHg and in fact both 
compounds are found to have similar cohesive energies. 
While BaMnHg may not be the most favorable target due 
to the stability of BaH2 as a competing phase, our results 



suggest that synthesis of (MnHg) 2 ~ salts may be possible 
and should be attempted. We note that the chemistry of 
Mn VII is different than that of Re VII. For example, 
KMnC>4 is highly oxidizing in aqueous solution, while 
KRe04 is not. Thus, while it may be possible to syn- 
thesize (MnHg) 2 ~ salts, the solvents and synthesis path 
needed may be different from that used in the (ReHg) 2 " 
salts. Light element salts of (MnHg) 2 ~ would have very 
high hydrogen contents, in excess of 10 weight % in the 
cases of Li, Mg or Be. 
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